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1. Introduction 

In a previous article [1], preliminary calculations of phase shifts were presented for 
the elastic scattering of low energy positrons by molecular hydrogen. The calculations 
used a generalized form of the Kohn variational method |2] and were carried out for 
the lowest partial wave of Eg symmetry, which has been shown [3] to be the only 
significant contributor to scattering processes for incident positron energies below 2 eV. 
The Kohn calculations were the first for the {e~^ — H2) system to treat the potential 
term corresponding to the H2 molecule explicitly. This introduced a complication [1] 
that had been avoided in earlier applications [31 HI [5] of the Kohn method by the use of 
the method of models [6]. 

Our implementation of the Kohn method has recently been made considerably 
more accurate. This is due most notably to the inclusion in the trial wave function of 
a greater number of terms describing short-range correlations between the molecular 
electrons and the positron, as well as to the use of improved methods for eliminating 
numerical anomalies due to so-called Schwartz singularities [71 [8] . We intend to publish 
details of our analysis of Schwartz singularities in a forthcoming article. In this paper, 
we will discuss results that have arisen during our efforts to improve the quality of the 
Kohn calculations by introducing a more flexible description of leptonic correlations 
than has previously been used. 

Our earlier calculations [1] included in the description of the electron-positron 
correlation a number of terms that were linear in the electron-positron distance. The 
importance of such terms was first demonstrated by Hylleraas [9] and their inclusion 
greatly increases the speed at which the description of the electron-positron interaction 
converges. As explained by Armour [10], this is due to the role of linear terms in 
the electron-positron distance in satisfying the Kato cusp condition [11]. The use in the 
earlier Kohn calculations of Hylleraas-type functions in the electron-positron coordinate 
was seen to increase significantly both the variational estimate, t]^, of the phase shift 
and the estimate of Zes, the positron annihilation parameter. The value of Zes measures 
the correlation between the positron and the target molecule and can be regarded as 
the effective number of molecular electrons available for annihilation with the positron. 
Methods for determining theoretical values of Z^s in the context of our Kohn calculations 
have already been outlined [1]. 

Discrepancies remained between the results of the earlier calculations and available 
experimental data, and our efforts to improve the theoretical model are ongoing. We 
have recently extended our Kohn calculations to assess for the first time the contribution 
to ?7v and Zes made by Hylleraas-type functions that are linear in the interelectronic 
distance. We have found that, when Hylleraas-type functions of this form are included 
in the description of the short-range leptonic correlations, the results of the Kohn 
calculations can be unreliable unless the wave function used to represent the target 
H2 molecule is itself made more accurate by the inclusion of Hylleraas-type functions in 
the interelectronic distance. Our observations are similar to those made by Van Reeth 
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and Humberston [12] in their Kohn calculations of s-wave phase shifts for positron 
scattering by atomic helium. Close-coupling calculations of electron scattering using 
inexact target states [131 E] have also shown that a balance must be maintained between 
the complexity of the treatment of target and scattering function electronic correlation. 

Our findings are important as the calculations that we will show to be unreliable 
can involve target states that account for nearly 97% of the correlation energy [15] and, 
coincident ally, give agreement with experimental results. When efforts are made to 
make the calculation more accurate, this agreement is no longer observed and the effect 
of including Hylleraas-type correlation functions in the interelectronic distance is found 
to be small. To cope with these difficulties, we present a numerical technique to analyze 
the variational results of calculations using inexact target states. This technique allows 
us to distinguish between reliable calculations and calculations which are unreliable due 
to an imbalance of target and scattering electronic correlation. 

2. Theory 

2.1. The complex Kohn variational method 

A description of the Kohn variational method is given elsewhere [2]. Calculations specific 
to the elastic scattering of positrons by molecular hydrogen have been carried out by 
Armour and coworkers [31 HI [5] and, more recently, by Cooper and Armour [1] . Only 
the essential aspects of the method will be repeated here. 

The objective of the calculation is to obtain an accurate approximation to the exact 
scattering wave function, from which it is straightforward to calculate variational 
estimates, rj^, of the phase shift, as well as estimates of Z^s. Approximate wave functions 
can be obtained by prescribing a trial wave function, ^'t, whose general form is known 
but which depends linearly on n unknown parameters. The basis of the Kohn method is 
the application of a stationary principle that allows optimal values of these parameters to 
be determined by solving a corresponding system of linear equations in the n unknowns. 

The complex Kohn method [16] is an extension of the original variational approach 
and allows the trial wave function, \E't, to be complex- valued. It has been shown [17] that 
the use of complex trial functions can mitigate the effects of anomalous results due to 
Schwartz singularities that arise when the linear system of Kohn equations is numerically 
ill-conditioned. Our own investigations of Schwartz singularities, to be discussed in a 
future article, have confirmed that the results presented here are free of this type of 
anomalous behaviour. 

In our calculations on (e"*" — H2) scattering, we have found it convenient to fix the 
nuclei at the equilibrium internuclear separation, R = 1.4 a.u. and to label the electrons 
as particles 1 and 2, taking the positron to be particle 3 and describing each particle by 
prolate spheroidal coordinates [H] {Xj, fij, j G {1,2,3}. We have used a complex 
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trial wave function of the form 

M 



'^t = {S + atT + poXo)^G + J2Pi^i^ (1) 



i=l 



where 



T = 5 + iC, (2) 

N 

S = y^^n[c{\s-l)] (3) 



and 



N 

C = cos [c (A3 - 1)] {1 - exp [-7 (A3 - 1)]}. (4) 

A3 — i 

The functions 5* and C are the same as those used in our previous Kohn calculations 
P and represent, respectively, the incident and scattered positrons asymptotically far 
from the target H2 molecule. The shielding parameter, 7, ensures that C is regular at the 
origin and, as before P, is taken to have the value 7 = 0.75. The constant c is defined 
to be c = kR/2, k being the magnitude of the positron momentum in atomic units. N 
is a normalization constant. The unknowns at and {po, ...,pm} are the complex-valued 
constants to be determined by the Kohn variational method. The function ipQ is an 
approximation to the ground state wave function of the unperturbed hydrogen molecule 
and is determined by the Rayleigh-Ritz variational method [15]. The general form of 
ipQ will be discussed in section 12.21 

The short-range correlation functions, Q = {xoi ■■■iXm}, allow for the description 
of direct electron-positron and electron-electron interactions, xo is the same correlation 
function used in our previous calculation [1] and was introduced first by Massey and 
Ridley [I9]. The general form of the remaining functions, for {1 < i < M), is 

X^ = N [Af A^>?7i*Si3(^.) + A^A^>*/i^2»S23(^0 

X A3>3'si2(^i) exp [-/3 (Ai + A2) - aAg] , (5) 

for prescribed basis states, {oj, 6j, q, di, ri, Si, 9i}, comprising non-negative integers. The 
interparticle functions, Spg{9i), have the form 

P12 = 1^12 (Oi = 1) 

Sl2m = { Afi2COs(0i-(/)2) {9, = 2) (6) 

1 (otherwise). 



and 



f'is = rTjs {9i = 3) 

s.sie,) = { M,3 cos{4>, - 03) {9, = 4) (7) 
I (otherwise), 
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for j G {1,2}, where rpg is the distance between leptons p and q. The inclusion of terms 
of the form Mpg cos(</)p — (f)q), where 



(A^-l)(l-/^p)(A?-l)(l-/i?) (8) 



is equivalent to considering terms in r^^. The choice of the nonlinear parameters, a and 
P, will be discussed in section 12.31 

We have carried out Kohn calculations using two different sets of correlation 
functions, which for convenience we shall denote by Q^^^ and Q^^l The set Q^^^ 
has M = 279 and, in addition to xo, contains three subsets of 87 basis functions 
corresponding to 6i = 0, 6i = 2 and 6i = 4, as well as 18 Hylleraas-type basis functions 
in the electron-positron coordinates, for which 9i = 3. ^l^^^ has the same general form 
as the set of 99 correlation functions used in our earlier calculation [1]. The set has 
M = 297 and is identical to but for the inclusion of a further 18 Hylleraas-type 
basis functions in the interelectronic coordinate, for which 6i = 1. Further details of the 
individual basis functions used are available from the corresponding author. 



2.2. The hydrogen molecule 

Although the Schrodinger equation for the hydrogen molecule cannot be solved exactly, 
very accurate numerical approximations to the exact solution can be obtained. A 
standard approach for determining approximate wave functions of bound states is the 
Rayleigh-Ritz variational method, used to great effect by James and Coohdge [20] and 
Kolos and Roothaan [21] in their pioneering calculations on the hydrogen molecule. 
Following these authors, we have taken the approximate wave function, i/jq, to have the 
form 

L 

v=l 

where 




X exp [-5 (Ai + A2)] , (10) 

for prescribed basis states, {m„, n^, j^,, /c^,, u;^,}, comprising non-negative integers. The 
function 312(00^) has the same definition as used in Optimal values of the unknown 
constants {c^} are determined in the Rayleigh-Ritz method by minimizing the energy 
expectation of ^g- 

The accuracy of ipc is typically measured in terms of the correlation energy of the 
molecule. This is the amount of the ground state energy, due to electron correlation, 
beyond that which is taken into account in a Hartree-Fock calculation [13]. The 
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percentage, Pq, of the correlation energy accounted for by an approximate target wave 
function with ground state energy expectation, -Ecaic? is 

TP TP ^ ' 

where E'ex is the exact nonrelativistic ground state energy in the Born-Oppenheimer 
approximation [22] and E'hf is the Hartree-Fock energy. 

We have carried out Kohn calculations using two different target wave functions, 
which for convenience we shall denote by ijj^^ and ^/'^f'*. The function tp'^'^ has L = 144, 
with a basis set comprising 72 terms having a;^ = and 72 terms having = 2. A value 
of (5 = 1.14 was chosen for to minimize its ground state energy expectation value, 
which accounted for 96.8% of the correlation energy of H2. ipl^^ has the same general 
form as the 31-term function used in our earlier calculations [1]. The function ipf^^ has 
a 145-term basis set of an identical form to that used for but for the inclusion 

of one Hylleraas-type term in pu for which Uy = 1. The value of 6 for ipQ remained 
fixed at 6 = 1.14, and the corresponding ground state energy accounted for 99.7% of 
the correlation energy of H2. Further details of the individual basis functions used are 
available from the corresponding author. The values of E'hf and E^x used to determine 
Pc were taken respectively from the calculations of Jensen [23] and Wolniewicz |24j . 

The important role played by Hylleraas-type functions in pi2 in describing electronic 
correlations in the hydrogen molecule has long been known [20]. However, until very 
recently it was not feasible for us to carry out Kohn calculations with target functions of 
this form, due to difficulties in evaluating the corresponding integrals found in the Kohn 
equations. However, we have successfully made modifications to the computational 
framework used in our calculations so that target functions containing Hylleraas-type 
terms in pi2 can now be admitted. Earlier changes to this framework had already 
been made during our previous calculations [1], where code designed originally for 
investigations of helium-antihydrogen scattering [25] was adapted so that it could be 
applied to (e^ — H2) scattering. Those initial modifications allowed for the evaluation 
of integrals containing terms in P13P23/P12 by using a triple Neumann expansion 



To carry out the calculations described here involving ipQ , it was necessary to extend 
these modifications to allow for the evaluation of integrals containing factors of the form 

Gu3 = F{Xs)^^, (12) 

Pl3 

and 

G213 = F(A3)^^, (13) 

P23 

where F (A3) can be either of the two open-channel functions (13]|^ representing the 
positron, or the function xo- 
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2.3. Optimization 

In contrast to variational calculations of bound states, there is no energy minimization 
principle associated with scattering wave functions. As a result, there is no absolutely 
rigorous method available to optimize the nonlinear parameters, a and (3, characterizing 
the rate of decay of the short-range correlation functions. Nevertheless, arguments for 
preferred choices of these parameters can be made. 

For atomic scattering, it has been shown [28] that, for a system where no bound 
state exists, the Kohn method gives an upper bound on the scattering length, a, where 



and hence a lower bound on the exact phase shift, r/, in the limit of zero positron energy. 
In the case of the Kohn variational method, obtaining bounds on scattering phase shifts 
is not generally possible at all incident energies considered, owing to the occurrence 
of the Schwartz anomalous behaviour at certain energies. However, an analysis of the 
method for potential scattering by Brownstein and McKinley [29] showed that, away 
from these energies, the phase shift will be bounded, provided the trial functions are, in 
some sense, sufficiently accurate. 

In the case of the solution of scattering systems using a close-coupling expansion 
[30] , it has been shown that bounds exist on scattering phase shifts or eigenphase sums 
provided that all open channels are included in the expansion and that the open channel 
target states are exact [311 [32]. If additional correlation functions are added to the 
expansion over target states, the bounds are still valid under certain conditions [311 [32]. 
the extra terms acting as an optical potential for channels not explicitly included. 

Kohn calculations using exact target states have been carried out by Humberston 
and are described, for example, in [26]. He found that, at low energies, the variational 
approximation to the phase shift tended to increase monotonically as the flexibility 
of the trial wave function was improved by the inclusion of a greater number of 
short-range correlation functions. He concluded that it was reasonable to expect the 
variational approximation to converge upwards to the exact phase shift with the use of 
an increasingly flexible trial wave function. In the method of models it is assumed that 
the target wave function used in the calculation is an exact solution of a model problem, 
so that the Kohn scattering parameters converge to the exact values for the model. 

In the case of inexact target states, there are no known rigorous bounds on 
scattering parameters for general close-coupling calculations or for the Kohn method. 
The experience in both low energy elastic positron scattering [12] using the Kohn method 
and close-coupling calculations of electron scattering (see, for example, [13], [H]) is that 
the monotonic increase in phase shifts or eigenphase sums may continue well above the 
physical values if the description of scattering electronic correlation is made noticeably 
more intricate than the description of target electronic correlation. This highlights the 
need for determining whether the target states are sufficiently accurate to give reliable 
values for scattering parameters in a given calculation. 




(14) 
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Returning to our own implementation of the Kohn method, we see that, assuming 
our calculations are reliable in the sense that we have described here, we may regard r]^ 
as an effective lower bound for rj. Under these circumstances, values of the nonlinear 
parameters, a and (3, can justifiably be chosen to maximize r]y. 

3. Results and discussion 

Our calculations were carried out for the lowest partial wave of S+ symmetry. A total of 
four different trial wave functions were used, corresponding to combinations of the two 
sets of correlation functions, fi*^^-' and and the two inexact target functions, ipQ^^ 
and V'g • ^ self-evident nomenclature, we will denote the four different trial wave 
functions by \E'^^''^'', '^t^'^\ "^[^'^^ and "^[^'^K Values of ?7v and Z^s were determined for 
each trial wave function, for positron momenta in the range k = 0.01 a.u. to = 0.4 
a.u., corresponding to a maximum positron energy of 2.18 eV. As noted in section [H 
higher partial waves become significant in scattering processes above this energy. 

Following our discussion in section [2731 we selected values of a and (3 approximately 
to maximize 77^. In principle, such maxima could be found, at least numerically, 
by straightforward iterative approaches. In practice, however, repeating our Kohn 
calculations for different values of a and /? is computationally very expensive; as we have 
discussed elsewhere [1], each iteration necessitates the evaluation of a large number of 
integrals that can be obtained only numerically via a triple Neumann expansion. Our 
analyses were therefore restricted to a relatively small set of candidate values for the 
nonlinear parameters, namely, a G {0.2, 0.3, . . . , 0.9, 1.0} and P G {0.2, 0.3, . . . , 1.4, 1.5}. 
Kohn calculations were performed for each of the 126 combinations of a and /3. Unless 
otherwise noted, all of the results presented here are from calculations carried out with 
a = 0.3 and (3 = 0.7, which we found to maximize rj^ for the trial wave function, ^E'l^''^^ 
at A; = 0.04. This value of k is approximately equal to the average momentum of a 
Maxwellian distribution of positrons at 297 K. It is convenient to consider the positron 
distribution at this temperature as it allows a direct comparison to be made of our 
results for Zcs with experimental data. 

3.1. Calculations involving tpQ^^ 

We consider first the two trial wave functions involving ^/'q^'', having Pq = 96.8. The 
dependence of rj^ and Zes on k for ^^^'^^ and \E'|^'^'' is shown in figures [H and O We have 
also included in these figures the values of rj^ and Z^s found in table 2(e) of the account of 
Kohn calculations made by Armour and Baker Those calculations used the method 
of models with a trial wave function containing M = 72 correlation functions, including 
eight Hylleraas-type terms in the electron-positron distance. 

The effect of including in '^f'^^ the Hylleraas-type correlation functions in pi2 is 
clear. The calculated values of both 77^ and Z^s for vl/^^'"^^ are significantly greater at low 
positron momenta than the corresponding values for \E'P''^\ The differences between the 
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results for and become smaller at higher positron momenta. There is broad 

agreement between our results for rj^ and those reported by Armour and Baker, although 
there is insufficient data available from those calculations to determine whether better 
agreement is observed for \[/|^'^'' or '^[^''^\ Estimates, a^, of the total scattering cross 
section could also be determined directly from values of rj^- However, there is a paucity 
of available experimental cross-section data at the very low positron momenta of most 
interest here, making any meaningful comparison with our results very difficult. 



There is good agreement between the estimates of Zcs for 



(1,A) ^{2,A) 



and the 



results of Armour and Baker, for k > 0.1. Again, a lack of available data from those 
earlier calculations prevents a comparison below k = 0.1, where the differences between 
the results for '^[^'^^ and vl/p'^-* are striking. We can, however, remark that the calculated 
value of Zes for vl/p'^'' at /c = 0.04 is 13.4, in reasonable agreement with the accepted 
experimental value of Z^s = 14.61 ± 0.14 at 297 K [33] . 

The extent of the influence on 77^ and Z^s of the Hylleraas-type functions in pi2 
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becomes even more pronounced if the values of the nonhnear parameters, a and /3, are 
varied. Figures [3] and H] illustrate the respective dependence of rj^ and Z^s on a and 
f3 at k = 0.04, for \['|^'^'' and '^[^'^\ The effects of the Hylleraas-type functions in pi2 
included in \|/[^'^^ are most obvious for a < 0.5 and become more dramatic as the value 
of a decreases. Indeed, the values of rj^ and Zgs for vEfl^'^-* shown in figures [3] and H] have 
not reached an obvious plateau with respect to further decreases in the value of a. It 
seems plausible that these values would continue to increase with decreasing a. 

In view of this, we think it necessary to examine the possibility that the observed 
effect is not genuine and instead arises from inaccuracies in the numerical evaluation 
of the integrals required to formulate the Kohn equations. These could occur because 
the short-range correlation functions become more diffuse as the value of a decreases, 
increasing the range of the configuration space of the positron over which the effects 
of the correlation functions are significant. To investigate this, we carried out a more 
detailed study of the Kohn calculations at {a,P) = (0.2,1.1) and {a,P) = (0.2,0.8), 
corresponding respectively to the maximum values of rj^ and Zes observed for vE^p'^-* in 
figures [3] and HI If the effects observed for "^[^'^^ are due to problems with convergence 
of integrals, increasing the range of integration in A should have a significant effect on 
the results for rj^ and Z^s- However, when we increased the range of our integration in 
A by 50%, the values of ?7v at (a,/?) = (0.2, 1.1) and Z^s at (a,/5) = (0.2,0.8) changed 
respectively by only 0.1% and 0.2% from the values shown for '^[^''^^ in figures [3] and 
m This is clear evidence that the effects we have described are not due to errors in the 
numerical integration. 

Figures [IHl] indicate that the apparent importance of the Hylleraas-type correlation 
functions in pi2 is a general feature of the calculation at low positron momenta. This 
is unexpected, since functions of this type do not address explicitly the key difficulty 
of describing correlations in terms of the electron-positron separation. Following our 
discussion in section 12.31 it is conceivable that the observed behaviour is a result of 
inaccuracies in the calculation due to the use of the inexact target wave function, V'g^^ 
despite its taking into account of 96.8% of the correlation energy. This claim is consistent 
with the findings of Van Reeth and Humberston [12]. In light of their conclusions, we 
investigated the sensitivity of our own Kohn calculations to changes in the accuracy 
of the target wave function. Basis functions were removed incrementally at random 
from creating a series of target wave functions of successively lower accuracies. 

After each removal, Kohn calculations were performed using the target wave function 
of reduced accuracy to determine values of 77^ and Z^s, each time for two trial wave 
functions having the sets of correlation functions ^l^^^ and A maximum of 70 basis 
functions were removed from the original set of 144 terms, at which point the target 
wave function accounted for 90.7% of the correlation energy of H2. The dependence of 
?7v and Zes on the accuracy, Pq, of each target wave function is shown in figures [5] and 
[6] respectively, for k = 0.04. 

The values of r]^ and Zcs tend to decrease with increasing Pq, for both ^l^^^ and 
fi*^^) . However, there is a distinct convergence of r/v and Z^g for the results corresponding 
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Figure 3. r/v {a, (3) at fc = 0.04, for l-*^'^^ (top) and f^'^^ (bottom). 



to above Pq ~ 96, which is not evident for the results corresponding to More 
generally, we have observed similar behaviour for analogous calculations performed at 
different values of k. We can reasonably conclude that Kohn calculations carried out 
using ipf^^ are reliable when Hylleraas-type terms in are omitted from the correlation 
functions, but are unreliable when functions of this type are included. For trial wave 
functions containing the set of correlation functions, we would expect any further 
increase in the accuracy of the target wave function to be accompanied by a noticeable 
drop in the values of r]^ and Z^fi- 

3.2. Calculations involving tpQ 

We consider now the two trial wave functions involving ipQ , having Pc = 99.7. The 
corresponding Kohn calculations could be performed only when the modifications to 
the computational framework, described in section 12.21 had been implemented. The 
dependence of rj^ and Z^s on k for and vE^p'^-* is shown in figures [7] and [H We have 

also reproduced in these figures the values of t]^ and Zes determined earlier for 
and 
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Figure 5. The dependence of r]v on Pc at fc = 0.04, for [xjfi^^^ and [+]ri(^^ Basis 
functions have been removed successively from i/jq''^ . 
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A number of comments should be made about our results. Firstly, there is excellent 
agreement in the values of rj^ and Z^s between calculations using ^1^'^^ and 
substantiating our earlier claim that i/jq"^ is a sufficiently accurate target wave function 
for the Kohn calculation using vl/l^'^-* to have converged and hence to be considered 
reliable. Secondly, there are significant differences in the results for \Ef{: ' and \l/t • 
The improvement in the accuracy of the target wave function has brought the results for 
\E^t^'^'' broadly into line with those for vEfl^'^-* and \E't^'^''. When the more accurate Kohn 
calculations using i}jq are carried out, therefore, the effect of including Hylleraas-type 
correlation functions in pi2 is small. 

It remains to be shown that ipQ is a sufficiently accurate target wave function for 
calculations involving ^!\' to be considered reliable. To do this, we again removed basis 
functions successively at random from the target wave function to reduce its accuracy. 
After each removal, Kohn calculations were performed using the target wave function 
of reduced accuracy to determine values of r]^ and Z^q, each time for two trial wave 
functions having the sets of correlation functions Vt^^^ and A maximum of 104 

basis functions were removed from the original set of 145 terms comprising V'g ' 
which point the target wave function accounted for 90.4% of the correlation energy of 
H2. The first 70 terms removed corresponded directly to the 70 terms removed earlier 
from V^ff'^ Thereafter, the remaining 34 terms were removed successively at random, 
with the condition that the Hylleraas-type term in pi2 was not removed. The dependence 
of ?7v and Zgfj on the accuracy, Pq, of each target wave function is shown in figures [9] 
and [To] respectively, for k = 0.04. 

The convergence of r]^ and Zcs with increasing Pc is not as obvious as it was in the 
previous calculations illustrated in figures [5] and [6l however, the slopes of the curves in 
figures [9] and fTOl are becoming noticeably fiat in the upper limit of Pc- We have already 
concluded that the Kohn calculations involving ^1^1^'^^ and hence those involving ^E't^'^^ 
are reliable. Further, the behaviour of the curves corresponding to and is 
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very similar in the limit of high Pq in both figures M and (TUl Inspection of these figures 
suggests that neither set of results would change significantly if the target wave function 
was again extended to account for most of the remaining 0.3% of the correlation energy. 
We may reasonably regard the calculations involving \E'{: ' as reliable. 

An interesting feature apparent from figures [9] and [10] is that the inclusion of the 
Hylleraas-type term in pu raises the threshold of convergence for trial wave functions 
containing In these figures, the values of r]^ and are still clearly declining at 
Pc = 96, at which value we have already concluded that Kohn calculations containing 
using target wave functions without Hylleraas-type terms, have converged. The 
origin of this effect is not clear and will remain a subject of our investigations. 

4. Concluding remarks 

We have demonstrated that the reliability of Kohn calculations for (e"*" — H2) scattering 
can depend upon the flexibility of the correlation functions used in the trial wave 
function, relative to the flexibility and the accuracy of the approximate wave function 
representing the target. This dependence is most prominent at very low positron 
momenta. Our flndings are similar to those reported by Van Reeth and Humberston for 
positron-helium scattering and highlight the need for rigorous testing of the accuracy of 
Kohn calculations whenever inexact target wave functions are used. 

We have implemented a numerical method to test the stability of any given 
calculation to variations in the accuracy of the approximate target state. This has 
allowed us to distinguish between reliable and unreliable results and thus compensate 
for the lack of explicit bounds on the scattering phase shifts. 

Having carried out the most accurate of our Kohn calculations, we have observed 
that the effect of including Hylleraas-type correlation functions in pi2 is to increase the 
calculated values of t]^ only slightly. The changes in the values of Z^s are also small. 
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so that there is still disagreement between our reported values and the established 
experimental result of Z^s = 14.61 ± 0.14 at 297 K [33]. This discrepancy is significant 
when compared to the results of other applications of the Kohn method for simpler 
systems. The best available calculations for atomic helium [31], for example, obtain a 
theoretical value of Z^s = 3.88±0.01 at 293 K. The experimental value is 3.94±0.02 |35j . 
Our intention is to address the problems in our Kohn calculations for H2 by improving 
the flexibility of the correlation functions still further to include, for example, terms 
linear in both pi2 and pjs, j G {1,2}, as well as terms describing virtual positronium 
formation. Consideration of virtual positronium has been shown [36] to enhance 
signiflcantly the calculated values of Zcs for positron scattering by atomic hydrogen 
near the positronium formation threshold. We hope that a similar increase in Z^s 
will be observed in our own calculations for molecular hydrogen if virtual positronium 
formation is taken into account. In any event, we will try to obtain converged results 
with as flexible a set of short-range correlation functions as possible. 
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